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1 Introduction 



In recent years many attempts have been made to understand the nucleon structure functions 
measured in lepton deep- inelastic scattering (DIS). Ahhough perturbative QCD is successful 
in describing the variation of structure functions with the squared momentum transfer, their 
magnitude and shape is governed by the non-perturbative physics of composite particles, and 
is, so far, not calculable directly from QCD. 

A variety of models has been invoked to describe nucleon structure functions. Bag model 
calculations for example, which are driven by the dynamics of quarks bound in a nucleon bag, 
quite successfully describe non-singlet unpolarized and polarized structure functions (see e.g. 

and references therein). However such calculations are not relativistically covariant. 

A covariant approach to nucleon structure functions is given by so called "spectator models" 

^, 0. Here the leading twist, non-singlet quark distributions are calculated from the process 
in which the target nucleon splits into a valence quark, which is scattered by the virtual photon, 
and a spectator system carrying baryon number 2/3. Furthermore the spectrum of spectator 
states is assumed to be saturated through single scalar and vector diquarks. Thus, the main 
ingredient of these models are covariant quark-diquark vertex functions. 

Until now vertex functions have been merely parameterized such that the measured quark 
distributions are reproduced, and no attempts have been made to connect them to some dy- 
namical models of the nucleon. In this work we construct the vertex functions from a model 
Lagrangian by solving the Bethe-Salpeter equation (BSE) for the quark-diquark system. How- 
ever, we do not aim at a detailed, quantitative description of nucleon structure functions in 
the present work. Rather we outline how to extract quark-diquark vertex functions from Eu- 
clidean solutions of the BSE. In this context several simplifications are made. We consider only 
scalar diquarks as spectators and restrict ourselves to the SU{2) flavor group. The inclusion 
of pseudo-vector diquarks and the generalization to SU (3) flavor are relatively straightforward 
extensions and will be left for future work. It should be mentioned that the quark-diquark 
Lagrangian used here does not account for quark confinement inside nucleons. However the use 
of a confining quark-diquark interaction should also be possible within the scheme that we use. 

As an important result of our work we find that the vertex function of the nucleon is highly 
relativistic even in the case of weak binding. Furthermore we observe that the nucleon structure 
function, Fi, is determined to a large extent by the relativistic kinematics of the quark-diquark 
system and is not very sensitive to its dynamics as long as the spectator system is treated as a 
single particle. 

The outline of the paper is as follows. In Sec.^ we introduce the spectator model for deep- 
inelastic scattering. Section |^ focuses on the scalar diquark model for the nucleon which yields 
the quark-diquark vertex function as a solution of a ladder BSE. In Sec.^ we present numerical 
results for the quark-diquark vertex function and the nucleon structure function, Fi. Finally 
we summarize and conclude in Sec.^. 
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2 Deep-Inelastic Lepton Scattering in the 
Spectator Model 

Inclusive deep-inelastic scattering of leptons from hadrons is described by the hadronic tensor 

W^^{P,q) = ^ fd'^e''^<{PmOrmP), (1) 

ZTT J 

where P and q are the four-momenta of the target and exchanged, virtual photon respectively, 
and J'^ is the hadronic electromagnetic current. In unpolarized scattering processes only the 
symmetric piece of W'^'^ = W^^ is probed. It can be expressed in terms of two structure 
functions Fi and F2, which depend on the Bjorken scaling variable, x = Q^/2P ■ q, and the 
squared momentum transfer, = — g^: 

VI'-(,.P)=(-,-+«^)fUx.Q^)+(p''-«"^) (p"-'!'^)^^^ (2) 

In the Bjorken limit (Q^, P -q 00; but finite x) in which we work throughout, both structure 
functions depend up to logarithmic corrections on x only, and are related via the Callan-Gross 
relation: F2 = 2xFi. 

The hadronic tensor, W^'^, is connected via the optical theorem to the amplitude T^'^ for 
virtual photon-nucleon forward Compton scattering: 

-lmT^"'{q,P) = W^"'{q,P). (3) 

In the Bjorken limit the interaction of the virtual photon with a valence quark from the target 
leads to a spectator system carrying diquark quantum numbers, i.e. baryon number 2/3 and 
spin or 1. In the spectator model it is assumed that the spectrum of spectator states can 
be saturated through a single scalar and pseudo- vector diquark |^, 0. In the following we 
will restrict ourselves to contributions from scalar diquarks only. The generalization to include 
a pseudo-vector diquark contribution is left for future work. The corresponding Compton 
amplitude is (Fig. [|): 

Tr(9,^) = {l + ^)Jj0i:HP)mp-mk)YS{k+q)YS{k) (4) 

xD{P - k)T{k, P ~ k)u{P), 

where the flavor matrix has to be evaluated in the nucleon iso-spin space. The integration 
runs over the quark momentum k. The Dirac spinor of the spin averaged nucleon target 
with momentum P is denoted by u{P). Furthermore S{k) = l/{mq — — ie) and D{k) = 
l/(m|) — k'^ — ie) are the propagators of the quark and diquark respectively, while F is the 
quark-diquark vertex function. To obtain the hadronic tensor, the scattered quark and the 
diquark spectator have to be put on mass-shell according to Eq.@: 

S{k + q) inSiml- {k + qY){mg + lit + i), 

D{P-k) i7T6{ml-{P-ky). (5) 
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The vertex function F for the target, which in our approach is a positive energy, spin-1/2 
composite state of a quark and a scalar diquark, is given by two independent Dirac structures: 



r(fc,p - fc)i(p_,p_^ = (/r(fc') + '^^f^k')) (6) 

where A.^^\P) = 1/2 + If/2M is the projector onto positive energy spin-1/2 states and M = 
\fP^ is the invariant mass of the nucleon target. Note that according to the on-shell condition 
in Eq.(ll) the scalar functions, f^J^^ will depend on k"^ only. 

From Eqs.(§ - ^ we then obtain for the valence quark contribution to the structure function 



6 2 /at 167r3 i-oo m2-A;2 



1 



x^(l-.+ (-^tflr^.)ggl! (7) 



+ 4 




— A;2 



M J V M J ml-k^ 



The upper limit of the fc^-jntegral is denoted by: 

kL. = ^(M'--^], (8) 




Note that k'^^^ —>■ —oo for x ^ 1. This implies that for any regular vertex function _p™^- — > 
for X 1 and thus the structure function automatically has the correct support. 

Since the spectator model of the nucleon is valence-quark dominated, the structure function 
F™^' in Eq.(|^ is identified with the leading twist part of Fi at some typical low momentum 
scale, /i^ ^ 1 GeV^. The physical structure function at large ^ /i^ is then to be generated 
via Q^-evolution. 

It should be mentioned that the Compton amplitude in Eq.(P and also the expression 
for the structure function in (J^) contain poles from the quark propagators attached to the 
quark-diquark vertex functions. From Eq. (P) it follows that these poles do not contribute when 
M < mn + mq. This condition is automatically satisfied if the nucleon is considered as a bound 
state of the quark and diquark, as done in the following. 

In the next section we shall determine the vertex function F, or equivalently f"'^ and 
from Eq.(|3) as solutions of a ladder BSE. 
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3 Scalar Diquark Model for the Nucleon 



We now determine the vertex function as the solution of a BSE for a quark-diquark system. 
We start from the following Lagrangian: 

C= Ipa {i^ - rriq) Ipa + 9^0*9^ 0a " ml)(f)l(f)a (9) 

where we have explicitly indicated SU{3) color indices but have omitted flavor indices. We 
restrict ourselves to flavor SU{2), where T2 is the symmetric generator which acts on the iso- 
doublet quark field, ip, with mass rriq. The charged scalar field, 0, represents the fiavor-singlet 
scalar diquark carrying an invariant mass m^- Similar Lagrangians have been used recently 
to describe some static properties of the nucleon, such as its mass and electromagnetic charge 
(see e.g. §, 0, i). 

The nucleon with 4-momentum P and spin S is described by the bound state Bethe-Salpeter 
(BS) vertex function F: 

F.T. (0 \TMx)<f),{y)\P, S) = 6abS{k)D{P - k)iT{k, P - k)u{P, S). (10) 

Here u{P, S) is the nucleon Dirac spinor and F.T. stands for the Fourier transformation. [] 
(Again we have omitted SU{2) flavor indices.) 

We will now discuss the integral equation for the vertex function F in the framework of the 
ladder approximation. 



3.1 Ladder BSE 

For the following discussion of the integral equation for the vertex function F we write the 
quark momentum as g + rjgP and the diquark momentum as —q + r/D-P- The weight factors rjg 
and rjD are arbitrary constants between and 1 and satisfy rjg + rju = 1. Within the ladder 
approximation the BSE for the vertex function of a positive energy, spin-1/2 model nucleon 
can be written as (see Fig.(|^)): 

F(g, P)u{P, S) = 9'J j0^^S{-k-q-{vg-VD)P)S{v,P+k)D{r^DP-k)T{k, P)u{P, S), (11) 

where the flavor and color factors have already been worked out. The scattering kernel is given 
by a M-channel quark exchange according to the interaction Lagrangian in Eq.(^. 

Since we are only interested in positive energy solutions, we may write the vertex function 

as: 

r(g,P)= (^a/i(g,P) + 6/2(g,P) + ^/2(g,P))AW(P). (12) 

The arguments of the scalar functions fa{q,P) are actually and P ■ q, but we use this 
shorthand notation for brevity. With a and b we denote as yet unspecified scalar functions 
of and P ■ q which will be chosen later for convenience. (The definition of Eq.(^ 
corresponds to a specific choice of a and b.) 

1 We use the normalization {P'\P) = 2P°{2'k)^6(^\P' - P) and Y.s "(^' S)u{P, S) = \/p2 + If^M + If. 
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3.2 Wick Rotation 



After multiplying the BSE in ([TID with appropriate projectors (which depend on a and 6), we 
obtain a pair of coupled integral equations for the scalar functions P) and /2(g', P)'- 

f^{q,P)=g^J^^ D,i~q-k~iv,-Vn)P) (13) 

xD^ir]^P + k) Dniv.P - k) K^p{q, k, P) fp{k, P), 

where Dq{p) = l/{mg — — ie) and Dd{p) = l/(m|, — — ie) are the denominators of the 
quark and diquark propagators, respectively. The indices a and (3 stand for the independent 
Dirac structures of the vertex function F, i.e. in the scalar-diquark model they run from 1 
to 2 according to (|T2|). Consequently the function Ka/3{q, k, P) is a 2 x 2 matrix, where its 
explicit form depends on the definition of the scalar functions fa{q,P)- We use a form factor 
for the quark-diquark coupling which weakens the short range interaction between the quark 
and the diquark and ensures the existence of solutions with a positive norm. For simplicity, we 
use a M-channel form factor which can be conveniently absorbed into the denominator of the 
exchanged quark propagator as follows: 

D,ip) ^ D,ip) ^ D,ip) ^\ . . (14) 

As a next step let us analyze the singularities of the integrand in Eq. (p!3D . For this purpose we 
choose the nucleon rest frame where P^ = P^^ = (M, ) and put the weight constants r]q and 
riD to the classical values: 

_ mq _l-ri 
Vq = = , (15) 

Vd = = ■ (16) 

Here we have introduced the asymmetry parameter rj = , such that the invariant quark 

and diquark mass is given by rriq = m(l — 77) and mu = m(l + rj) respectively, where rh = 
{niq + m£))/2. In the complex ko plane Dq{riqP + k) and Dij{riqP — k) will be singular for: 

k' =-r]qM±Eq{k)T^e, (17) 
k^ =riDM ±EDik)Tie, (18) 



where Eq{k) = ym'^ + k and E^lk) = y mj^ + k . The cuts lie in the second and forth 
quadrant of the complex fco-plane. However for a bound state, < M < mq + niD, a gap occurs 
between these two cuts which includes the imaginary ko axis. 

Next, consider the singularities of the exchanged quark propagator: 

k^ = + r]M± Em, (g + k) t ^e, (19) 
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where ErnXk) = y mf + k and rrii = m,, A for i = 1,2, respectively. The sum of the second 
and third term at the RHS of Eq. (p!9D is bound by: 

T]M + EmAq + k) >{mD-m,)^^^ + mi, (20) 
r]M-EmXq + k) <{mD-m,)j^^-m,. (21) 

The diquark should be considered as a bound state of two quarks which implies rriD < 2mg. 
Together with rrii > niq, namely setting the form factor mass A larger than rUg, we have rriu — 
mq+rrii > and mD — rrig — mi < 0. Consequently we find from Eqs.(^,|2T]) rjM +Emi{q+k) > 
and rjM — E„ii (q + k) < for any momenta q and k. Therefore, if — g° + r^M — Emi (q + k) > or 
—q^ + rjM + Emiiq + k) < a so-called "displaced" pole will occur in the first or third quadrant, 
respectively. In other words, the displaced-pole-free condition is: 

7]M - E.^Xq + k)<q'' <r]M + EmXq + k), (22) 

for any k. Since k is an integration variable, EmXq + k) will adopt its minimum value, rriq, at 
k = —q for i = 1. The above condition therefore simplifies to: 

(g° - riMf < mj. (23) 

If g° is Wick rotated to pure imaginary values, i.e. q^^ q^^ = {iq'^, q) with real G (— oo, oo), 
the displaced poles will move to the second and forth quadrant. Then, after also rotating the 
momentum k^ ^ k^ = (ik^, k), we obtain the Euclidean vertex functions fp{k,P^^^) from the 
Wick rotated BSE: 

/„(g,pW)=^?2/0 f),i-q-k-iv,-r^n)P^'^) (24) 

Z},(r/,P(°) + ~k) Do{v,P^'^ - k) K^p{q, k, P^) fp{k, p(°)), 

where d^kE = dk'^d^k. 

If we are in a kinematic situation where no displaced poles occur, i.e. Eq.(|2^) is fulfilled, we 
may obtain the Minkowski space vertex function /□-(g, P) from the Euclidean solution through: 

^(g,pW) = ^?2/0 ^^(_5_^_(^^_^^)p{0)) (25) 

xD,(r/,pW + k) Do{v,P^''> - k) K^^iq, k, P^) fp{k, P^). 

It should be emphasized that for a given Euclidean solution fi3(k,P^^^), Eq. (^5|) is not an 
integral equation but merely an algebraic relation between fa{q, P^^-*) and fjs^k, P^^^). If however 
displaced poles occur, i.e. Eg . (^3]) is not fulfilled, one needs to add contributions from the 
displaced poles to the RHS of Eq. (P3|) . This will lead to an inhomogeneous integral equation 
for the function faiq, P^^^), where the inhomogeneous term is determined by the Euclidean 
solution fp{k,P^'^^). 
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Since the Euclidean solutions /^(g, P^°^) are functions of = —q% = —((54)^ + 
q . p(o) = iq^M for a fixed M, it is convenient to introduce 4-dimensional polar coordinates: 





= qE COS aq 




= \q\q', 




= qE sin ttq 



Here < < tt and the 3- dimensional unit vector is parameterized by the usual polar and 
azimuthal angles = (sin6'g cos(/)g, sin 0g sin 0q, cos6'g). In the following we therefore consider 
fa{q,P^'^^) as a function of qE and cosa^. Furthermore it is often convenient (and traditional) 
to factor out the coupling constant together with a factor (47r)^, and define the "eigenvalue" 
= {g/Any. Then the BSE in (p^ ) is solved as an eigenvalue problem for a fixed bound 
state mass M. 



3.3 0(4) Expansion 

In the following we will define the scalar functions fa{q, P) for positive energy (M > 0) bound 
states via: 

P-q 



P) 



h{q.P) + 



A(+)(P), 



(27) 



i.e., we now make a specific choice for the scalar functions a and h in Eq.(|T^). In the rest frame 
of this model nucleon, this leads to: 



<:l',r (q^ p''') = nq, p^'^HP^'\ s) = ( ij^'^pll^ ) X.. 



(28) 



Here we have explicitly used the Dirac representation. The Pauli matrices a act on the two 
component spinor xs: where the spin label S* = ibl is the eigenvalue of a^: cr^xs = Sxs- In 
terms of the 0(3) spinor harmonics yf^^ [^], 



y'oi'iq) 



Xs 



and yl^^ {q) = -q. ay^^^^ (g) , 



(29) 



we have: 



=1/2+ 

J3=S/2 



{q,P^'^) 



fiiq,P^'^)yol"iq) \ 

-^Sf2{q,P^'^)yl^\q)J- 



(30) 



From Eq.(|30D we observe that /i and /2 correspond to the upper and lower components of the 
model nucleon, respectively. 

After the Wick rotation, as discussed in the previous subsection, the scalar functions fa 
become functions of qE and cosa^. Therefore we can expand them in terms of Gegenbauer 
polynomials Cl{z) 



fa{qE, cos aq) = J2 ^"/a('?i?)Cn(cos 



On 



(31) 



n=0 
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X 



We have introduced the phase to ensure that the coefficient functions are reaL The 
integral measure in 0(4) polar coordinates is: 

W)'=w?k Jo ^"'^ 2^^ y ^''^ 

Multiplying the BSE in Eq.(piD with the Gegenbauer polynomial C^(cosa;g) and integrating 
over the hyper-angle a^, reduces the BSE to an integral equation for the 0(4) radial functions 

-Pn. 
J a ■ 

\{M)f2{qE) = Y.ll / dkElCl^{qE.kE) f}^{kE). (33) 

13=1 m=0 •'^ 

Here A(M) is the eigenvalue which corresponds to a fixed bound state mass M. Furthermore 
note that the integral kernel 

2 /"^ 2 /"^ 

lC^'g{qE,kE) = {-lY i"^ - dag siii^ ag Cl{cos ag) - / rfofc sin^ C^(cosafc) 

^ TT Jo 71 Jo 

I- J dn^ k% Dgi-q -k-ir]g- r]D)P^'^) K^p{q, k P^) (34) 

X D,(r/,P(0) + fc) Do{v,P^^^-k) 

is real, so that we can restrict ourselves to real 0(4) radial functions 

To close this section we shall introduce normalized 0(4) radial functions. Since the scalar 
functions fi{q,P) and /2(g, P) correspond to the upper and lower components of the model 
nucleon respectively, one may expect that /2(g, P) becomes neghgible when the quark-diquark 
system forms a weakly bound state. Thus one can use the relative magnitude of the two scalar 
functions, f2{q, P)//i(g, P), as a measure of relativistic contributions to the model nucleon. To 
compare the magnitude of the Wick rotated scalar functions fi{q, P) and f2{q, P) we introduce 
normalized 0(4) radial functions. Recall the 0(4) spherical spinor harmonics [|1^, 0: 

Z„,,^(a,^,0)= n{smayClt\icosa)yUOA). (35) 

The integers n and / denote the 0(4) angular momentum and the ordinary 0(3) orbital angular 
momentum respectively. The half-integer quantum numbers j and m stand for the usual 
0(3) total angular momentum and the magnetic quantum number. We rewrite the Wick 

rotated solution {q,P^^^) in terms of the spinor harmonics Znjim{ce,9,(l)) and define 

the normalized 0(4) radial functions P„(g_B) and Gn{qE) as: 

The extra factor \/2 n is introduced for convenience. The normalized 0(4) radial functions P„ 
and Gn are then linear combinations of the /^: 

FniqE) = fUQE), (37) 

/ \ r~, ( f2~^{(lE) . f2^^{.(lE)\ /oo\ 

Gn{qE) = -^\/-(- + 2)^^ + ^^ • (38) 
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Equivalently we can express the Wick rotated scalar functions as: 



fi{qE,cosa) = ^ C^(cosa), (39) 



n=0 



/2(g£;,cosa) = . (7^(9^) Cn-i(cosa). (40) 

n=i Jn(n + 2) 



3.4 Euclidean Solutions 

In this section we present our results for the integral equation in (|33|). For simplicity we 
considered the quark and diquark mass to be equal: rriq = mj:, = rh. In this case the kernel 
JC"^^ can be evaluated analytically in a simple manner, since for rj = the denominator of the 
propagator for the exchanged quark does not depend on the nucleon momentum P. We fixed 
the scale of the system by setting the mass m to unity. The "mass" parameter in the form 
factor was fixed at A = 2m. 

We solved Eq. (^3]) as follows. First we terminated the infinite sum in Eq. (|3lD at some fixed 
value nmax- Then the kernel in Eq. (|3^ ) for the truncated system becomes a finite matrix 
with dimension (2 x ^max)^- Its elements are functions of qe and ks- Next we discretized the 
Euclidean momenta and performed the integration over numerically together with some 
initially assumed radial functions f^. In this way new radial functions and an "eigenvalue" A 
associated with them were generated. The value of A was determined by imposing the normal- 
ization condition on such that resultant valence quark distribution is properly normalized. 
We then used the generated radial functions as an input and repeated the above procedure 
until the radial functions and A converged. 

Note that our normalization differs from the commonly used one [|^, since we are going 



to apply the vertex function only to processes with the diquark as a spectator, i.e., we do 
not consider the coupling of the virtual photon directly to the diquark. The ordinary choice 
of normalization would not lead to an integer charge for the model nucleon in the spectator 
approximation. We therefore normalize the valence quark distribution itself. 

Regarding Eq.(^) as an eigenvalue equation, we found the "eigenvalue" A (coupling con- 
stant) as a function of M^, varying the latter over the range 0.85 < < 1.99 m^. The 
eigenvalue A was stable, i.e. independent of the number of grid points and the maximum value 
of kE- Furthermore, for a weakly bound state, M > 1.6 m, the solutions were independent 
of the choice of the starting functions and the iteration converged in typically 10 ~ 25 cycles. 
However, for a strongly bound state, M < 1.4 m, we found that the choices of the starting func- 
tions were crucial for convergence. A possible reason of this instability for a strongly bound 
system is due to the fact that we did not use the 0(4) eigenfunctions F„ and G„ in numerical 
calculations but the functions defined in Eq. (0). Since strongly bound systems, M ~ 0, 
are approximately 0(4) symmetric, a truncated set of functions may be an inappropriate 
basis for numerical studies of the BSE. 

We found that the eigenvalue A converges quite rapidly when rimax, the upper limit of the 
0(4) angular momentum, is increased. This stability of our solution with respect to n^^^ is 
independent of M. We observe that contributions to the eigenvalue A from failE) with n > 4 
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are negligible. This dominance of the lowest 0(4) radial functions has also been observed in 
the scalar-scalar-ladder model and utilized as an approximation for solving the BSE in a 
generalized fermion-fermion-ladder approach [p!5[] . 



To compare the magnitude of the two scalar functions P) and f2{q, P) we show in Fig.^ 
the normalized 0(4) radial functions F„ and As the dependence of A on rimax suggests, 
radial functions with 0(4) angular momenta n > 4 are quite small compared to the lower 
ones. Together with the fast convergence of A, this observation justifies the truncation of 
Eg . (|33|) at n = nmax- Note that even for very weakly bound systems (~ 0.5% binding energy) 
the magnitude of the "lower-component" f2{q,P) remains comparable to that of the "upper- 
component" fi{q,P). This suggests that the spin structure of relativistic bound states is 
non-trivial, even for weakly bound systems. So-called "non-relativistic" approximations, in 
which one neglects the non-leading components of the vertex function {f2{q, P) in our model), 
are therefore only valid for extremely weak binding, 2m M only. 

3.5 Analytic Continuation 

In the previous section we obtained the quark-diquark vertex function in Euclidean space. Its 
application to deep-inelastic scattering, as discussed in Sec.^, demands an analytic continuation 
to Minkowski space. Here the scalar functions fa, which determine the quark-diquark vertex 
function through Eq.(0), will depend on the Minkowski space momenta and P ■ q. 

Recall that our Euclidean solution is based on the expansion of the scalar functions fa in 
terms of Gegenbauer polynomials in Eq.(^). This expansion was defined in Sec.^]^ for real 
hyper-angles ag, with —1 < cosa^ < 1. Consequently the infinite sum over the 0(4) angular 
momenta n in (|3l|) is absolutely convergent for pure imaginary energies q'^. Now we would like 
to analytically continue g° to physical, real values. The Euclidean hyper-angle aq is defined in 
Euclidean space such that: 

COSCtg = . — (41) 

In Minkowski space cos ag is then purely imaginary (cos ag = —i (f j \/~q^) for space-like g, and 
real (cosa^ = —q^/\fq^) if q is time-like. Note that the angular momentum sum (|3T|) converges 
even for complex values of cosa^ as long as | cosa^l < 1. Then an analytic continuation of fa 
to Minkowski space is possible. In terms of the Lorentz invariant scalars q^ and P -q we obtain: 

z = cos ag = -sgn(g2) J^^ . (42) 

Then the convergence condition for the sum over the 0(4) angular momenta in Eq. (|3l|) reads: 

{P ■ q^ < \q^\. (43) 

Even if Eq.(^) is satisfied the radial functions themselves may contain singularities 
which prevent us from performing the analytic continuation by numerical methods. However, 
note that the Euclidean solutions for f^ are regular everywhere on the imaginary g°-axis. 
Consequently the RHS of the "half Wick rotated" equation (El) contains no singularities if 
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the displaced-pole-free condition in Eg .(1231) is met. Therefore, in Minkowski space, the radial 
functions fa are regular everywhere in the momentum region where the displaced-pole-free 
condition (|23|) and the convergence condition (^) are satisfied. Here the analytic continuation 
to Minkowski space is straightforward. Recall the normalized 0(4) radial functions F„ and Gn 
from Eqs. ([37|j38|) , which are linear combinations of f^- Writing them as Fn{qE) = Qe ^niQE]) 
and Gn{qE) = Qe Gn{qE), we find for the scalar functions fai^f, P ■ l) from Eqs. (|39|j40| ): 



2\ 



/l(?^^•?) = E^^ (V^^M^j Cliz), (44) 



n=i ^Jn{n + 2) 



/2(g^P■g) = -E / mJ. (V^^M^j Cl_,{z). (45) 



Note that the Gegenbauer polynomials (O^-i) together with the square root factors (V?^ M^)" 
((Vg2 A-f^)"^^) are n-th ((n — l)-th) order polynomials of g^, M^, and P ■ q and contain there- 
fore no \/q^ M'^ factors. Since in the kinematic region under consideration, /i(g^,P ■ q) and 
/2(g^,P ■ q) are regular, it is possible to extrapolate Fn{—q^) and Gn{—q^) numerically from 
space-like q^ to time-like q^ as necessary. 

Finally we are interested in the quark-diquark vertex function as it appears in the handbag 
diagram for deep-inelastic scattering. Therefore we need the functions fa for on-shell diquarks 
only. The squared relative momentum q^ and the Lorentz scalar P ■ q are then no longer 
independent but related by: 



F ■ q = 



2m 



D 



\mq + mD) ^ 



(46) 



Then fi and /2 from Eq. (0) and (|45|) are functions of the squared relative momentum q^ only. 

In Sec.0 the parameterization (^) for the Dirac matrix structure of the vertex function was 
more convenient to use. The corresponding functions f^ which enter the nucleon structure 
function in Eq.(^ are given by : 

fr{k') = \f2{q\P-q). (48) 

Here the arguments q^ and F ■ q, of the scalar functions fa on the RHS should be understood 
as functions of fc^ through Eq. (^Bj) , together with the relation: 

= (k^ - m, 7710]]. (49) 

As already mentioned, the procedure just described yields radial functions /°" in Minkowski 
space only in the kinematic region where the conditions Eqs.(^,^) are met. These are ful- 
filled for weakly bound states < {niD + at moderate values of On the other 
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hand, the nucleoli structure function in (|^ at small and moderate values of x is dominated 
by contributions from small quark momenta < m^. Consequently, the Minkowski space 
vertex function obtained in the kinematic region specified by the displaced-pole-free condition 
( p3D and the convergence condition (|^) determines the valence quark distribution of a weakly 
bound nucleon at small and moderate x. 

In the case of strong binding, <^ {mg+mDY , or at large x the nucleon structure function 
is dominated by contributions from large space-like k"^. Here the above analytic continuation 



to Minkowski space is not possible and the sum over the 0(4) angular momenta in Eg . (|3TD 
should be evaluated first. In principle this is possible through the Watson-Sommerfeld method 
jl6], |l^ where the leading power behavior of fa{q^,P ■ q) for asymptotic P ■ q can be 
deduced by solving the BSE at complex 0(4) angular momenta , or by assuming conformal 
invariance of the amplitude and using the operator product expansion technique as outlined in 
ref.fg. 

However, the use of the operator product expansion is questionable here, since in the quark- 
diquark model which is being used we have introduced a form factor for the quark-diquark 
coupling and our model does not correspond to an asymptotically free theory. Existence of the 
form factor also makes the analysis of complex 0(4) angular momenta complicated. Therefore 
a simpler approach is used. It can be shown from a general analysis that BS vertex functions 
which satisfy a ladder BSE are regular for space-like fc^, when one of the constituent particles 
is on mass shell. Furthermore, from the numerical solution studied in the previous section, we 
found that the magnitude of the 0(4) partial wave contributions to the function decreases 
reasonably fast for large 0(4) angular momenta n, except at very large k"^. We therefore use 
the expansion formulae (^4D and (^5]) with an upper limit on n < ^max to evaluate /°° defined 
by Eqs.( |47| , ^8|) as an approximation. Nevertheless, this application of BS vertex functions to 
deep-inelastic scattering emphasizes the need to solve Bethe-Salpeter equations in Minkowski 
space from the very beginning, as has been done recently for scalar theories without derivative 
coupling pO . 



4 Numerical Results 

In this section we present results for the valence contribution to the nucleon structure function, 
Fi, from Eq.(0), based on the numerical solutions discussed above. First we show in Fig. Ijthe 
physical, on-shell scalar functions for a bound state mass M = 1.8 m. The maximal 0(4) 
angular momentum is fixed at n^ax = 4. Figure § demonstrates that the magnitude of /{"^ and 
/2'^ is quite similar, even for a weakly bound quark-diquark system. Furthermore we find that 
for weakly bound states (M > 1.8 m) the dependence of on rimax is negligible in the region 
of moderate, space-like — fc^ ^ 5 m^. However for larger space-like values of fc^ the convergence 
of the 0(4) expansion in Eq. (|3lD decreases for any M^, and numerical results for fixed rimax 
become less accurate. 

In Fig. ^ the structure function, F™'', is shown for various values of using rimax = 
4. The distributions are normalized to unity. One observes that for weakly bound systems 
(M = 1.99 m) the valence quark distribution peaks around x ~ 1/2. On the other hand. 
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the distribution becomes flat if binding is strong (M = 1.2m). This behavior turns out to 
be mainly of kinematic origin. To see this, remember that F™^' is given by an integral (c.f., 
Eq (0)) over the squared quark momentum, k^, bounded by fc^^^^ = a;(M^ — m|,/(l — x)). The 
latter has a maximum at x = 1 — niD/M. Therefore the peak of the valence distribution for 
weakly bound systems occurs at x ~ 1/2 for = mq. For a more realistic choice tud ~ 2mq 
the valence distribution would peak at x ~ 1/3. The more strongly the system is bound, the 
less fc^ax varies with x. This leads to a broad distribution in the case of strong binding. Thus 
the global shape of F™'' is determined to a large extent by relativistic kinematics. 

To investigate the role of the relativistic spin structure of the vertex function we discuss 
the contribution of the "relativistic" component, /I", to the nucleon structure function F™^'. 
Figure ^ shows that the contribution from is negligible for a very weakly bound quark- 
diquark state (M = 1.99m). Here the "non-relativistic" , leading component, /°°, determines 
the structure function. However, even for moderate binding the situation is different. In Fig.^ 
one observes that the contribution from the "relativistic" component is quite significant for 
M = 1.8 m. Nevertheless, the characteristic x dependence, i.e. the peak of the structure 
function at x ~ 1/2, is still due to the "non-relativistic" component. 

5 Summary 

The aim of this work was to outline a scheme whereby structure functions can be obtained from 
a relativistic description of a model nucleon as a quark-diquark bound state. For this purpose 
we solved the BSE for the nucleon starting from a simple quark-diquark Lagrangian. From 
the Euclidean solutions of the BSE we extracted the physical quark-diquark vertex functions. 
These were applied to the spectator model for DIS, and the valence quark contribution to the 
structure function Fi was calculated. 

Although the quark-diquark Lagrangian used here is certainly not realistic, and the cor- 
responding BSE was solved by applying several simplifications, some interesting and useful 
observations were made. We found that the spin structure of the nucleon, seen as a rela- 
tivistic quark-diquark bound state, is non-trivial, except in the case of very weak binding. 
Correspondingly, the valence quark contribution to the structure function is governed by the 
"non-relativistic" component of the nucleon vertex function only for a very weakly bound 
state. Furthermore, we observed that the shape of the unpolarized valence quark distribution 
is mainly determined by relativistic kinematics and does not depend on details of the quark- 
diquark dynamics. However at large quark momenta difficulties in the analytic continuation of 
the Euclidean solution for the vertex function to Minkowski space emphasize the need to treat 
Bethe-Salpeter equations in Minkowski space from the very beginning. 

Acknowledgments 

This work was supported in part by the the Australian Research Council, BMBF and the 
Scientific Research grant #1491 of Japan Ministry of Education, Science and Culture. 



14 



References 

[I] R.L. Jaffe, Phys. Rev. Dll, 1953 (1975). 

[2] A.W. Schreiber, A.W. Thomas and J.T. Londergan, Phys. Lett. B237, 120 (1989); Phys. 
Rev. D42, 2226 (1990); A.I. Signal and A.W. Thomas, Phys. Rev. D40, 2832 (1989). 

[3] H. Meyer and RJ. Mulders, Nucl. Phys. A528, 589 (1991). 

[4] W. Melnitchouk, A.W. Schreiber and A.W. Thomas, Phys. Rev. D49, 1183 (1994), Phys. 
Lett. B335, 11 (1994). 

[5] W. Melnitchouk, G. Piller and A.W. Thomas, Phys. Lett. B346, 165 (1995). 

[6] A. Buck, R. Alkofer and H. Reinhardt, Phys. Lett. B286, 29 (1992). 

[7] S. Huang and J. Tjon, Phys. Rev. C49, 1702 (1994). 

[8] H. Asami, N. Ishii, W. Bentz and K. Yazaki, Phys. Rev. C51, 3388 (1995). 

[9] A.R. Edmonds, Angular momentum in quantum mechanics, (Princeton Univ. Press, 1957). 

[10] A. Erdelyi, et. al, Higher Transcendental Functions (McGraw-Hill Book Company Inc., 
New York, 1953). 

[II] K. Rothe, Phys. Rev. 170, 1548 (1968). 
[12] K. Ladanyi, Nucl. Phys. B57, 221 (1973). 

[13] N. Nakanishi, Suppl. Prog. Theor. Phys. 43, 1 (1969). 

[14] E. zur Linden and H. Mitter, Nuovo Cim. 61B, 389 (1969). 

[15] P. Jain and H.J. Munczek, Phys. Rev. D48, 5403 (1993). 

[16] M.L. Goldberger and K.M. Watson, Collision Theory (Wiley, New York, 1964). 

[17] G. Domokos, Phys. Rev. 159, 1387 (1967). 

[18] M.L. Goldberger, D.E. Soper and A.H. Guth, Phys. Rev. D14, 1117 (1976). 

[19] M.L. Goldberger, D.E. Soper and A.H. Guth, Phys. Rev. D14, 2633 (1976). 

[20] K. Kusaka and A. G. Williams, Phys. Rev. D51, 7026 (1995). 



15 



Figure Captions 

Fig. ^ The diquark spectator contribution to the virtual forward Compton amplitude in the 
Bjorken limit. 

Fig. 0: The Bethe-Salpeter equation for a quark-diquark system in the ladder approximation. 

Fig. ^: The normalized 0(4) radial functions F„ and G„ from ( |36|) for M = 1.8 m as functions 
of the Euclidean momentum qe- 

Fig. ^: The on-shell scalar functions (solid) and —f^'^ (dashed) as a function of the quark 
momentum, k"^, for M = 1.8 m and nmax = 4. 

Fig. ^ The valence quark distribution, F™'', from Eq. (|^) for different binding for the model 
proton. The solid, dashed, and dot-dashed lines show the results for weak (M = 1.99 m), 
moderate (M = 1.8m), and strong (M = 1.2 m) binding respectively. 

Fig. 1^: Contribution of the "relativistic component", /j"^, to the structure function F™^' from 
Eq. (^) for a weakly bound model proton (M = 1.99 m). The solid line denotes the total valence 
distribution, F™^'. The dashed line shows the contributions to F™^', which are proportional to 

son 
J2 ■ 

Fig. 1^: As in Fig.|^, but for moderate binding for the model proton (M = 1.8 m). 
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